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We obtain information on the QED photon amplitudes at high orders in perturbation theory starting from known results on the QED effective 
Lagrangian in a constant electric field. A closed-form all-order result for the weak field limit of the imaginary part of this Lagrangian has 
been given years ago by Affleck, Alvarez and Manton (for scalar QED) and by Lebedev and Ritus (for spinor QED). We discuss the evidence 
for its correctness, and conjecture an analogous formula for the case of a self-dual field. From this extension we then obtain, using Borel 
analysis, the leading asymptotic growth for large A'^ of the maximally helicity violating component of the / - loop A*' - photon amplitude in 
the low energy limit. The result leads us to conjecture that the perturbation series converges for the on-shell renormalized QED N - photon 
amplitudes in the quenched approximation. 
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1. Introduction: The perturbation series in 
quantum electrodynamics 

In quantum electrodynamics, physical quantities are usually 
computed as a perturbative series in powers of the fine struc- 
ture constant a — j;^: 

F{e^) ^ Co + C2e^ + cte^ + . . . (1.1) 

In the early years of quantum field theory it was hoped that 
this type of series would turn out to be convergent in gen- 
eral, or at least for sufficiently small values of the coupling 
constant(s). However, in 1952 Dyson [1] argued on phys- 
ical grounds that QED perturbation theory should be diver- 
gent. Dyson argues: "Suppose, if possible, that the series 
(1.1) converges for some positive value of e^; this implies 
that f (e^) is an analytic function of e at e = 0. Then for 
sufficiently small values of e, F(— e^) will also be a well- 
behaved analytic function with a convergent power-series ex- 
pansion" [1]. He then argues physically that this cannot be 
the case, since for < the QED vacuum will be unstable 
due to a runaway production of e+e~ pairs which coalesce 
into like-charge groups. This argument does not prove diver- 
gence or convergence, but it gives a first hint of the connec- 



tion between perturbative divergence and instability in quan- 
tum field theory [2-5]. 

To substantiate this divergence mathematically in the case 
of QED turns out to be exceedingly difficult [2]. Explicit all- 
order calculations are, of course, next to impossible, so that 
the natural thing to try is to prove the divergence by estab- 
lishing lower bounds on the contributions of the individual 
Feynman diagrams representing the coefficients of the pertur- 
bation series. However this is still very hard to do for QED 
diagrams since here the integrands have no definite sign. On 
the other hand, in scalar field theories the integrands have a 
definite sign (in the Euclidean) so that here lower bounds can 
be established using inequalities such as 

At about the same time that Dyson's paper appeared. Hurst 
[6] and ThiiTing [7] (see also Petermann [8]) used this lower 
bound strategy to show by explicit calculation that perturba- 
tion theory diverges in scalar Xcj)'^ theory for any value of the 
coupling constant A. Hurst's proof is quite straightforward; 
its essential steps are (i) the use of (1.2) in the parametric rep- 
resentation to find lower bounds for arbitrary diagrams (ii) a 
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proof that the number of distinct Feynman diagrams at n* 
loop order grows like {'^)\n\ (iii) the fact that there are no 
sign cancellations between graphs. This led Hurst to con- 
jecture that the analog statement holds true for other renor- 
malizable quantum field theories and in particular for QED. 
However, QED was already well-established experimentally 
in 1952. To account for this fact, Hurst postulated that the 
QED perturbation series, while not convergent, is an asymp- 
totic series and thus still makes sense numerically. 

Today the fact that the perturbation series is asymptotic 
rather than convergent is not only well-established for QED, 
but believed to be a generic property of nontrivial quantum 
field theories (see [2-5, 9, 10]). Convergence of the series 
can be expected only in trivial cases where higher-order ra- 
diative corrections are absent altogether (usually on account 
of some symmetry). 

In the absence of convergence, the reconstruction of the 
exact physical quantities from their perturbation series must 
be attempted using summation methods. For this there are 
various possibilities, and we mention here only the one which 
has been most widely used in QFT, Borel summation [9, 10]. 
For a factoriaUy divergent series 

oo 

F(a)^^c„a"+i (1.3) 
one defines the Borel transform as 

°° i" 

B(t) = Vc„- (1.4) 
■^-^ n! 

n=0 

Assuming that B{t) has no singularities on the positive real 
axis and does not increase too rapidly at infinity, one can also 
define the Borel integral 

POO 

F{a) = / dte-^B{t) (1.5) 
Jo 

F is the Borel sum of the original series F. F is asymptotic to 
F by construction, although the physical quantity represented 
by the series F might still differ from F by nonperturbative 
terms. Even when the Borel transform has singularities it re- 
mains a useful concept, since these singularities contain de- 
tailed information on the divergence structure of the theory. 
In many cases they can be traced either to instantons (related 
to turmeling between vacua) or Euclidean bounces (related 
to vacuum decay) or to renormalons. Those latter typically 
arise in renormalizable theories and are related to large or 
small loop momentum behaviour. Diagramatically, they can 
be analyzed in terms of "infinitely long" chains of 'bubble' 
diagrams (for a review, see [11]). 

Despite of the many insights which have been gained 
along these Unes, a point which remains poorly understood is 



the influence of gauge cancellations on the divergence struc- 
ture of a gauge theory. Generally, in gauge theory individual 
diagrams do not give gauge invariant results; gauge invari- 
ance is recovered only after summing over certain classes of 
diagrams. In QED a textbook example is provided by one- 
loop photon-photon scattering, where a gauge invariant (as 
well as UV finite) result is obtained only after performing a 
sum over the six inequivalent orderings of the external mo- 
menta for the basic diagram in fig. 1. 




Fig. 1: Sum of one loop photon scattering diagrams. 

This recovery of gauge invariance generally implies can- 
cellations in the sum over gauge related diagrams, and leads 
one to expect that the coefficients of the perturbation series 
for such amplitudes should come out smaller in magnitude 
than predicted by a naive combinatorial analysis. However, 
it is presently not known what implications this has for the 
large order behaviour of amplitudes in either QED or other 
gauge theories. In 1977 Cvitanovic [12] performed a de- 
tailed diagramatic analysis of the effect of gauge cancella- 
tions in a previous calculation of the sixth order contribution 
to the electron magnetic moment [13]. Based on this ex- 
ample, he suggested that in QED their effect is sufficiently 
strong to modify estimates based on the basic factorial growth 
in the number of Feynman diagrams; more realistic estimates 
might be obtained by counting gauge invariant classes of dia- 
grams instead. The asymptotic growth of the number of these 
classes is, however, less than factorial. For the case at hand, 
the electron g — 2,he conjectured that its perturbation series 
even converges in the 'quenched' approximation, i.e. for the 
contribution represented by diagrams not involving electron 
loops. 

In the present paper, we present an analogous conjecture 
for the N - photon amplitudes, in scalar and spinor QED. 
Here the 'quenched' approximation corresponds to taking 
only the diagrams involving just one electron loop, which is 
also the 0{Nf) part of the amplitude in QED with Nf flavors. 
Our conjecture is that this part of the amplitude converges 
in perturbation theory when renormaUzed on-shell. We use 
known results on the QED effective Lagrangian to obtain in- 
formation on the large order behaviour of this amplitude for 
the special case of 'all +' polarizations, in the limit of low 
photon energies and large photon number N. 
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2. The Euler-Heisenberg Lagrangian and the N 
- photon amplitudes: general case 

Let us start with some basic facts about the QED effective 
Lagrangian in a constant external field F^,^. At one loop in 
spinor QED, this is the well-known Euler-Heisenberg ('EH') 
Lagrangian, obtained 1936 by Heisenberg and Euler [14] in 
form of the following integral: 



jC^^'^. (F) = 



-— r 



dT 



{eaT){ebT) 



tanh(eaT) tan(e6T) 



(2.1) 



Here T is the proper-time of the loop fermion, m its mass, 
and a, b are the two Maxwell invariants, related to E, B by 
a? - b"^ = B'^ - E'^ , a6 = E • B. A similar representation 
exists for scalar QED [15,16]: 



Then 



T^^\kr,ei;...;kN,SN\=C^^\iFtot^x) \f,. 



(2.5) 



On the right hand side it is understood that, after expand- 
ing out C^^^ («-Ftotai) to the A''th order in the field, only those 
terms are retained which involve each i^i , . . . , Ffq linearly. 
Using this expansion together with a convenient choice of 
polarizations, one can obtain a closed-form expression for 
the low energy limit of the one-loop on-shell photon ampli- 
tudes, valid for arbitrary N and all polarization components 
[18, 19]. We will be concerned here only with the 'all -|-' or 
'maximally helicity violating' component of the amplitude, 
defined by choosing the polarization vectors for all photons 
as positive helicity eigenstates (the 'all -|-' and 'all — ' compo- 
nents are related by a parity transformation). This component 
becomes particularly simple: 



,(1) 

■'seal 



167r2 



dT 

r3 ' 



{eaT)iehT) 
sinli(ear) sin(e6r) 

(2.2) 



+ ±(a2 _ 62)r2 - 1 



By standard QFT (see, e.g., [17]) these effective La- 
grangians contain the information on the one-loop photon S- 
matrix in the limit where all photon energies can be neglected 
compared to the electron mass, Wj << m. Thus diagra- 
matically /^gp],^ (F) is equivalent to the sum of the Feynman 
graphs in fig. 2 at zero momentum. 




+ 



Fig. 2: Diagrams equivalent to the EH Lagrangian. 

To obtain the one-loop N - photon amplitude from the 
Lagrangian, introduce for each photon leg the field strength 
tensor 



Fr = fcr< - ^r^r (2.3) 

with £j, ki the photon polarization and momentum. Define 



N 



-Ftotal — Fi 



(2.4) 



Kl:dki,et;---;kN,s+] = -2 

clpin('^) = - 



(2e) 



N 



2n{2n - 2) 



c^l(f)Xiv 



(2.6) 



Here the Bk are the Bernoulli numbers, and xn is a kinemat- 
ical factor which in spinor helicity notation (see, e.g., [20]) 
can be written as 



(ill 
2f 



{[12]2[34]2...[(iV-l)iV]2 



permutations I 
(2.7) 



([ij] = {kf\k~)). Its precise form is not essential for our pur- 
pose. Here and in the following it should be understood that 
we deal with the low energy limits of the amplitudes only. 



3. Two loop: the self-dual case 

The two-loop correction C^'^\F) to the Euler-Heisenberg 
Lagrangian, involving an additional photon exchange in the 
loop, has been considered by various authors starting with 
V.I. Ritus in 1975 [21-24]. However, for a general con- 
stant field no representation is known for C'^'^'^ [F) which 
would be exphcit enough to derive closed all-iV formulas for 
the two-loop photon amplitudes. Things simplify dramati- 
cally, though, if one specializes the field to a (Euclidean) self 
dual ('SD') field, obeying 



F„, 



-£ rF'^'^ 



(3.1) 
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The square of F becomes then proportional to the identity 
matrix, F"^ = — /^l, and we can use the number / to 
parametrize F. f can be either real or imaginary. The real 
/ case corresponds to a field which in Minkowski space has 
a real magnetic component and an imaginary electric com- 
ponent. In Euclidean space (with 3:4 = it) the field strength 
matrix can be Lorentz transformed to the form 



0/00 



-/ 

(3.2) 

/ 



0-/0 



Thus it describes a "doubling up" of the magnetic case, and 
indeed the EH Lagrangian for such a field turns out to have 
the same qualitative properties as the EH Lagrangian for a 
purely magnetic field. The self-dual case with a real / will 
therefore be called 'magnetic' in the following. Similarly, for 
a self dual field with a purely imaginary / the properties of 
the Lagrangian are similar to the electric EH Lagrangian and 
we call this case 'electric'. 

In this self dual case all integrals turn out to be elemen- 
tary, leading to the following simple closed-form expression 
for the two-loop effective Lagrangian [18], 



■^ipin('^) 



-2a 



m 1 



[3€'(«)-^'(«)] 



(3.3) 



Here k ; 



and 



(3.4) 



c^^^. (n) 



(47r)2m' 
1 r 2n - 3 



(27r)2 \ 2n - 2 



2n-2 



+3E 



B2k B 



^ 2k (2n-2fc) 



(3.5) 



(the kinematic factor xn is the same as in (2.6)). 



4. Higher loop orders: Im C 

Beyond two loops, it gets rather hard to obtain information on 
the EH Lagrangian by a direct calculation. Nevertheless, the 
work of Ritus and Lebedev quoted above [21, 22] leads to an 
all-order prediction for the imaginary part of this Lagrangian. 
As is well-known, the EH Lagrangian will have an imaginary 
part whenever the field is not purely magnetic. Physically 
it represents the possibility of vacuum pair creation by the 
electric field component. At one loop for the case of a purely 
electric field, Schwinger [16] found for it the following rep- 
resentation. 



4 oo 



fe=i 



A:2 



Trfc 



(4.1) 



with /3 = Here the term with index k describes the co- 
herent production of k pairs by the field. This representation 
is manifestly nonperturbative in the field and coupling, mak- 
ing it clear that the imaginary part cannot be seen unless one 
sums up the whole set of diagrams in fig. 2. 

At the two-loop level, the imaginary part of £{E) was 
first studied by Lebedev and Ritus [22] who found the fol- 
lowing generaUzation of Schwinger's formula (4.1): 



^/3^ ^aT:Kk{l3) exp 



fe=i 



irk 



(4.2) 



where i(^{x) = T'{x)/r{x) is the digamma function. Since 
self dual fields are helicity eigenstates [25] the effective ac- 
tion for such fields carries precisely the information on the 
'all -I-' (or 'all — ') photon ampUtudes. Thus at the two-loop 
level we are still able to write down a closed-form a\l-N ex- 
pression for this particular polarization choice: 



Here a = |^ is the fine-structure constant. The exponen- 
tial factors are the same as in (4.1), however at two loop the 
kth exponential comes with a prefactor K). which is it- 
self a function of the field. These prefactor functions are not 
known explicitly, but [22] were able to show that they have 
small (3 expansions of the following form: 



4 



ci =0, 



.^ + 1 + 0(7^) 
fe-i 

i 

Cfe 



J_ V 

^Vk^^ ^i{k-iy 



k>2 



(4.3) 



In the following we will concentrate on the weak field limit 
« 1 of Im£. In this hmit the k - series in (4.1), 
(4.2) are dominated by the fc = 1 term, and from (4.3) 
lim/3_>o^^i(/3) = 1. Thus one has 



Im£ 



(4.4) 



Now [22] find that the an term in (4.4) has a natural in- 
terpetation in the pair creation picture. Naively, to turn real 
a virtual electron-positron must separate out along the field 
direction a distance ry such that the energy gained from the 
field makes up for the rest masses. This is 



rii = 



2m 



(4.5) 



However, this does not take into account the Coulomb attrac- 
tion; a pair getting "born" at a separation of r|| comes with a 
negative binding energy of — a/ry. Lebedev and Ritus argue 
that, since this energy reduces the amount of energy which 
has to be taken out of the field, it is equivalent to a lowering 
of the rest mass. Thus m, the physical (on-shell renormal- 
ized) vacuum mass, should be replaced by an effective field- 
dependent mass m*{E), 



m^{E) = m — 



a eE 
2 



(4.6) 



with corrections expected to be of higher order in /3 (see 
[26,27] for a detailed discussion). Then replacing m by 
mt,{E) in the one-loop k = 1 exponential and expanding 
in a one finds that 



exp 



— TT 



mm 

eE 



1 + aw + 0{a^) 



exp 



— TT 



m 



(4.7) 



Thus if this interpretation is correct then the 1 -|- qitt in the 
two-loop formula (4.4) is the truncation of an exponential 
series, and the missing higher order terms must show up in 
the higher loop corrections to Im£spin(£^). Moreover, in the 
weak field limit these should be the only missing terms, lead- 
ing to a remarkably simple all-order prediction: 



1=1 



87r3 



■ exp 



TT 



(4.8) 



In the weak field limit the above considerations are spin- 
independent [27] so that the same formula (4.8) with just 
an additional factor of 1/2 on the right hand side applies 
to scalar QED. For this case, moreover, there exists a com- 
pletely different and more direct derivation of the same for- 
mula due to Affleck et al. [28]. (These derivations were 
independent of one another: apparently Affleck et al. were 
not aware of Ritus's work, and Ritus and Lebedev were not 
aware of Affleck et al.'s work.) It uses the following 'world- 
line path integral' representation of the quenched effective 
action in scalar QED [29], 



p(quenched) ^^-j 



dT 



-m-'T 



J Vx{t) e-^[^(^)l 



(4.9) 



Here the path integral 'Dx{t) runs over the space of all 
embeddings x{t) of the particle trajectory into (Euclidean) 
spacetime, with periodic boundary conditions in proper-time 
x{T) = x{0). The path integral action S'[a;(r)] has three 
parts. 



So 

Se 
Si 



S — So + Se + Si 



(4.10) 



= le 



ie I x''A^{x{t)) 
Jo 

-^lo'^'Io 



dro 



x{ti) ■ x(t2) 
{x{T,)-x{T2)f 



(4.11) 



describes the interaction of the scalar with the external 
field. Si the exchange of internal photons in the scalar loop. 
In [28] this representation was apphed to the constant electric 
field case, and \mCsca\{E) was calculated using a stationary 
path approximation both for the T integral and the path inte- 
gral / 'Dx{t). The stationary trajectory turns out to be a cir- 
cle with a field dependent radius. Evaluation of the worldhne 
action S'[a;(r)] on this trajectory yielded exactly the exponent 
in (4.8), and the second variation determinant gave precisely 
the same prefactor. Affleck et al. then argued that the sta- 
tionary point approximation becomes exact in the weak field 
Umit up to renormalization effects. 
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Formulas analogous to (4.1), (4.2) exist also for the self 
dual case [30]. In the 'magnetic' case (real /) the EH La- 
grangian is real, while in the 'electric' case (imaginary /) it 
has an imaginary part. At one and two loops, this imaginary 
part has the following Schwinger type expansion. 



-27rkK 



m 1 



27r/cfc 



(4.12) 



Here the functions Kk{K,) are the SD analogues of the 
Lebedev-Ritus functions Kk{(i). In contrast to the Kk{(i), 
the K}^{k) are known explicitly as power series in 1/k [31]: 



Kk{K) 



1 o °° 

27r/t •TT ^ 



m=l 



(4.13) 



In the weak field limit / ^ or k ^ oo one finds the same 
relation between the one and two loop contributions to Im£ 
as in (4.4), 



a7rlm4p| 



n(»«) 



(4.14) 



Although neither the Lebedev-Ritus arguments nor the ap- 
proach of Affleck et al. carry over to the self dual case in an 
obvious way, we take this as evidence that the simple expo- 
nentiation (4.8) applies also to the self dual case. Thus we 
conjecture that 



m 1 



(47r)2 K 



exp 



-2ttk + aTT 



(4.15) 

5. Re>C ^ Im>C via Borel dispersion relations 

We would now like to convert the result (4.8) for Im£ into in- 
formation about Re£ and the photon amphtudes. The appro- 
priate dispersion relation involves the same Borel technique 
which we described in the introduction for the loop expan- 
sion, now apphed to the weak field expansion. Let us first 
demonstrate it at the one loop level [32]. In the purely mag- 
netic case, it is easy to obtain from (2.1) the following closed 
formula for the coefficients of the weak field expansion of 



C 



(1) 

spin' 



(B) 

spin V / 



n=0 



(5.1) 



where g 



\^m2 ) 



and 



22"B. 



'2n+4 



(2n-F4)(2n-|-3)(2n-F2) 



(5.2) 



Using properties of the BernouUi numbers one can then show 
that asymptoticaUy the behave as 



,(1) 



(-i)"r(2n + 2) 



1 



1 



22n+4 ' 32n+4 



(5.3) 



If one replaces On ^ with just the leading term in this expan- 
sion, the series (5.1) turns into an alternating series which is 
divergent but Borel summable, i.e. the Borel integral (1.5) 
is well-defined. The weak field expansion of -Cgpin for the 
purely electric case differs from (5.1) only by an additional 
factor (-1)": 



(E) 



2m^ 



E(-i)' 



(5.4) 



n=0 



where now g = [j^j ■ In this case the same leading ap- 
proximation produces a non-alternating series for which the 
Borel integral diverges. Nevertheless, the imaginary part of 
this integral is well-defined through a dispersion relation (us- 
ing the discontinuity across the cut along the negative g axis) 
and yields just the fc = 1 term in the Schwinger expansion 
of Im£gp|j^, eq. (4.1). Applying the same procedure to the 
whole asymptotic expansion (5.3) term by term one repro- 
duces the complete Schwinger expansion. 

At the two loon level, no closed formula is known for the 
corresponding coefficients a„ ' . Still, in [33] their leading 
and subleading asymptotic growth could be determined by 
calculating the first fourteen coefficients and using them for 
a numerical fit. This yielded 



.(2) 



(-l)"r(2n + 2)^ 



0-44 \ 



87r3 



r2n 



(5.5) 



The leading term differs only by a factor a-K from the cor- 
responding one loop term, and thus when used in the Borel 
dispersion relation gives the air term in (4.4). From the sub- 
leading term one obtains the second term in the small /3 ex- 
pansion (4.3) for Ki, 



(5.6) 



ifi(/3) = l-(0.44)W-/3-fO(/3) 




Thus the Borel technique allows one to map the asymptotic 
large n expansion of the weak field expansion coefficients a„ 
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of the magnetic Lagrangian to the weak field expansion of 
the imaginary part of the electric Lagrangian. 

In the self dual case, where closed formulas are available 
for the real and imaginary parts of the EH Lagrangian also at 
two loops, this correspondence can be checked in even more 
detail. In [30] we calculated, up to exponentially suppressed 
terms, the large n expansion of the two loop weak field ex- 
pansion coefficients c'^\ (n) for scalar QED (the c'^\ {n) dif- 

(2) 

fer from the cl^^{^{n) given in (3.5) only in the coefficient of 
the second term in braces, which in the scalar QED case is 
3/2 instead of 3). Using a Borel dispersion relation as above 
then yielded precisely the small / expansion of (the scalar 
QED equivalent of) Ki{k). 



1=1 



1=2 



1=3 





6. Large N behaviour of the i-loop N- photon 
amplitudes 

Assuming that this Borel dispersion relation remains valid at 
higher loop orders, we can apply it in reverse to the exponen- 
tiation formula (4.15) and obtain all-loop information on the 
coefficients of the weak field expansion: 



lim„ 



(6.1) 



By (2.5) this translates into a statement on the 'all +' ampli- 
tudes in the hmit of large photon number N, 



1=4 





lim 



JV— »oo 



r«Jfci,£+;...;fc..,4] (^-1)! 



(6.2) 



Surmning this relation over I we formally get an exponentia- 
tion formula for the complete amplitude in the large N limit: 



-i(total) 



lim 



[ki,s 



+. 

1 5 • ■ 



JV- 



+ 1 



(6.3) 



Assuming sufficient uniformity in I of the convergence of the 
ratio (6.2) for N ^ oo, one could now conclude that the 
amplitude must be analytic in a for some sufficiently large 
N. But analyticity of the complete amplitude is excluded 
by renormalons and other arguments. Therefore uniformity 
must fail, and it is easy to see how this comes about diagra- 
matically. Fig. 3 shows the Feynman diagrams for the N 
photon amphtude for the first four loop orders (the external 
photons are not displayed, as well as diagrams differing from 
those depicted only by moving the end points of photon lines 
along the fermion loops). At I loops, it involves diagrams 
with up to Z — 1 fermion loops. 



Fig. 3: Feyrmian diagrams up to four loop order. 

From the Affleck et al. approach to the derivation of the 
exponentiation formula (4.8) we know that the weak field 
limit oflmC{E) comes entirely from the 'quenched' part of 
the amplitude. Thus at fixed loop order I the total contribu- 
tion of all non-quenched diagrams must be subleading in A'' 
compared to the quenched ones. Since the number of non- 
quenched classes of diagrams grows with increasing loop 
number, clearly one would expect the convergence in N to 
slow down with increasing I. 

If, on the other hand, one stays inside the quenched class 
of diagrams from the beginning, there is no obvious mathe- 
matical reason to expect such a slowing down. Thus we be- 
lieve that, for the quenched amplitude, the perturbation series 
indeed converges for sufficiently large A'^, and in fact for all 
N > 4; since it would be very surprising to find a quaUtative 
change in the convergence properties of these amplitudes at 
some definite value of A^ (although such a phenomenon can- 
not be excluded altogether, of course). Similarly, considering 
the essentially topological character of gauge cancellations 
one would not assume the restriction to zero momenta or the 
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choice of polarizations to be essential. Thus we are led to the 
following generalization of Cvitanovic's conjecture: the per- 
turbation series converges for all on-shell renormalized QED 
amplitudes at leading order in Nf . 

It must be emphasized that on-shell renormalization is 
essential in all of the above. The finiteness of the hmits in 
(6.1) and (6.2) requires the weak field expansion coefficients 
Cgpjjj (n) to have the same leading asymptotic growth for all 
I, namely T{2n — 2) in the magnetic case and r(2n — 1) in 
the 'magnetic' self dual case. At the two loop level, one can 
show by explicit calculation (numerically for the magnetic 
case [33], analytically for the 'magnetic' self dual one [30]) 
that the use of any renormaUzed electron mass other than the 
on-shell mass leads to an asymptotic growth faster than at 
one loop. In fact, a simple recursion argument shows that, if 
mass renormaUzation is done generically, then the I loop co- 
efficients will grow asymptotically like r(2n + Z — 3) in the 
magnetic and hke r(2n -|- / — 2) in the 'magnetic' self dual 
case. 

7. Conclusions 

To summarize, we have used existing results on the imagi- 
nary part of the QED Euler-Heisenberg Lagrangian together 
with Borel analysis techniques and a number of explicit two 



loop checks to argue for convergence of the perturbation se- 
ries for the physically renormalized QED photon amplitudes, 
to leading order in Nf. To the best of our knowledge, this 
statement is not at variance with known arguments against 
convergence for the full amplitudes, since those generally 
rely on the presence of an infinite number of virtual particles. 
In fact, it presently remains an open possibility that conver- 
gence might hold to any finite order in Nf. 

It would be clearly desirable to further corroborate our 
conjecture by higher loop calculations. While an explicit 
three-loop calculation for the magnetic or electric EH La- 
grangian seems presently out of the question, for the self dual 
case this calculation might be technically feasible [34]. Since 
the three loop amplitude has already a non-quenched part, we 
can make for it the double prediction that (i) the weak field 
limit of the quenched contribution to lm£gpjj^(iK) should dis- 
play the (a7r)^/2 correction implicit in (4.15) and (ii) the 
non-quenched contribution should be suppressed in this limit. 
Moreover, one would hke to see the ratio c'^^\n) / d^\n) 
to converge not slower than c^'^\n) / c^^\n) if only the 
quenched contribution is used in calculating c^^^ (n). 
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